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Abstract. 

We discuss a new approach to reducing excited state contributions from two- and three-point correlation functions in lattice 
simulations. For the purposes of this talk, we focus on the A( 1232) resonance and discuss how this new method reduces excited 
state contamination from two-point functions and mention how this will be applied to three-point functions to extract hadronic 
form factors. 

PACS: 11.15.Ha,12.38.Gc,13.40.Gp,14.20.Gk 

The calculation of the electromagnetic form factors for mesons and baryons is crucial to understanding the structure 
of hadronic states in QCD. However, they are notoriously difficult to both measure experimentally as well as calculate 
theoretically due to the complications that arise from the strong interactions. In the case of the A(1232) resonance, 
the form factors themselves are not currently experimentally accessible, although two of the form factors in the static 
limit are known (the charge) or measured to some degree (the magnetic dipole moment). For the nucleon, experimental 
results do exist for the form factors as a function of the momentum transfer. This makes a calculation of the nucleon 
form factors both a check of methodology as well as of QCD, and allows us to be confident that lattice results for the 
A form factors are reasonable. 

The electromagnetic form factors of the A are encoded in the matrix element 

(A(p')\J^\A(p)) = u a {p')T aft „up{p) (1) 

where u a is a Rarita-Schwinger vector-spinor describing the external A, and = L^g/^g is the vector current. 
The Lorentz structure of F is given, for example, in [1], and has four form factors F^i^Q 2 ) that are functions 
of Q 2 — —{p 1 — p) 2 alone. These form factors give rise, in the limit Q 2 — )• 0, to the electric charge, magnetic dipole 
moment, electric quadrupole moment, and magnetic octupole moment of the A. 
Of these moments, the charge is of course known, and from the PDG [2], we have 

Ai A++ = (5.6±1.9)mjv, (2) 
Ma+ = (2.7±3.5)mjv, (3) 

where we have added all of the errors (including theoretical) in quadrature just to get an idea for how well these are 
determined experimentally. Thus, it is essential even for this simple quantity to have a well-determined lattice result. 
Some unquenched results were obtained using the form factor approach in [3] and using a background field technique 
[4], but there are several difficulties that arise in these different methods. For this work, we will focus on the difficulties 
with the form factor approach. 

These are determined by calculating the 3 -point correlator: 

C 3 P\t I JJ fl p il p f )=FT[(0\x(tf,x f )Mt,0)x(ti,Xi)\0)] , (4) 

where FT is the Fourier Transform of the correlator, and % is some appropriate interpolating operator for the A. In the 
large time limit f/ 3> t ^> 

C^M/.pi.p/) ^Z(p i ,p / ) e - £ /( f /-') e - £ --^)(A(p / )|^(0)|A^)> + --- • (5) 

Here we have schematically written this so that Z contains various overlap factors of the form (0|^|A) as well as 
other kinematic factors that are known. The dots denote contributions from excited states that are generally ignored. 
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FIGURE 1. On the left is the matrix element Ge for the A ++ evaluated at q = 0. While there is significant excited state 
contamination, there is a noticeable plateau that gives (up to renormalization) the charge of the A. On the right is the same matrix 
element for q 2 = 1 \(2z/a s L) 2 ~ 4 GeV 2 , where there is no trustworthy signal. 



The standard approach then is to note that the overlap factors and kinematics can be canceled by an appropriate ratio 
with 2-point correlators, and from these ratios we can extract the matrix element of interest. Note that additional 
complications arise in separating the A state and the N — it state when the pion mass is below 300 Me V. Currently we 
are not near this regime, so it is not an issue for the current analysis. 

The problem is that the contamination from excited states can be seen to be large, even for the simplest cases. Take 
here the A ++ £0 form factor, GEo{q 2 ), which is a linear combination of the F* form factors. In the limit q 2 — > 0, 
Ge(0) — -f 2, the electric charge of the A ++ in units of \e\. We show on the left of Fig. 1 the appropriate ratio of 3-to-2- 
point correlators to get this quantity using the Hadron Spectrum collaboration lattices of 2+1 -flavor anisotropic Clover 
lattices [5] (here with a volume of 16 3 x 128 and a pion mass of roug hly 390 MeV, and a' 1 w 5.5 GeV, a s /a, w 3.5). 
The source and sink are located att = 0,28, and we can see a plateau in the center where we could reliably extract the 
form factor, but there is significant contamination from excited states, seen from the deviation from the plateau. 

This problem is more pronounced at higher momenta, and we show the same form factor at q 2 — 1 1 (2n/a s L) 2 w 4 
GeV 2 on the right of Fig. 1. The location of a plateau is questionable at best for this plot. 

We can extract the form factor using this standard approach and one gets reasonable results. For the q 2 = point, 
we find the charge of the A ++ is precisely twice that of the proton, and we can extract the renormalization factor for 
the vector current (note it is not conserved here because we are using a local current) from Geo(0), and we obtain 

— = 1.05(1) (statistical errors only). 

Zy 

While this determination is trustworthy, extracting the higher momenta form factors is dangerous, due to the excited 
state contamination. So we would like to examine a new approach, which makes use of the Variational method to 
better extract states that contribute to a correlator. For now we will discuss this in the context of 2-point correlators for 
simplicity. 

For a given state, there are many interpolating operators that could be used to calculate a two-point function, and 
one could form a matrix of correlators 

Cij(t) = (0\Oi(t)dj(0)\0) . (6) 
From this, by solving the generalized eigenvalue problem (GEVP) 

C(f)x = A(?,?o)C(? )x, (7) 

one can show (see [6] and references therein) that the eigenvalues behave like 

l i (t,t ) ~e- m > ( f-^ + - . (8) 
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Using 2 operators, t0=4 




FIGURE 2. On the left we show the effective mass for a single A operator with noticeable excited state contamination before 
reaching a plateau. On the right is the effective masses for the two states extracted using the GEVP and two different local A 
operators. 



The reference time fo is empirically chosen so that at that time, all n states (for an n x n system) would contribute to 
the correlator; no more, no less. For our purposes, we find that the choice of fo has little effect on our results, primarily 
because we are using a very small basis of operators. 

If we restrict ourselves to local operators, there are two A interpolating fields we can use, and we can see that this 
has little effect to reduce contamination from excited states in the ground state, as shown by comparing the two plots 
in Fig. 2. On left is a single operator effective mass, and we can see the excited state contamination before the plateau. 
On the right is the case with two operators, and solving the GEVP. The ground state is unchanged, and although this 
allows us to perhaps extract an excited state, this is not what we are currently interested in. 

So we would like to find a way to kill off the excited state contribution in the ground state. This is where the 
Generalized Pencil-of-Function (GPoF) method comes in [7, 8]. In a quantum mechanical system, the important point 
is that if (9a (?) is an interpolating operator for the A, then so is 

Ol(t)^e m O A (t)e- m = A (t + x), (9) 

and this new operator is linearly independent from the original operator. 

So if we use A (t),0^(t) as our two operators, we can construct a correlator matrix using only a single correlator. 
This matrix is 



CO) = 



(O a (0<5a(0)> <O A (r)<5 A (0)> \ _ ( (O A (t)O A (0)) <O A (f + T)0 A (O)> 

(o A (0oj(0)> (oi(t)di(o)) J { (o A (t)d A (-T)) (o A (t + T)d A (-T)) 

C(t) C(t + z) 
C(f + r) C(t + 2x) 



(10) 



We can replicate this and use a set of operators A "(?) = A (t + nz), and make this correlator as large as we wish. It 
turns out that using multiple shifts does not give us any additional information, so we will ultimately use just a single 
shift and a T = 4. This would be something that is determined for each correlator one is interested in. Once this matrix 
is created then we follow through the same procedure as before with the GEVP. 

Note that when the correlator basis grows, there is more of a likelihood of zero singular values, so we perform an 
SVD cut on the correlator matrix to exclude states that have significantly smaller singular values than the largest by 
some cutoff, here chosen to be k 10~ 3 . 

In Fig. 3 we show the same operator(s) as before, but now using the GPoF method (again, z = 4). On the left is one 
operator with one shift and the right uses both local operators and a single shift. We see the effective mass comes to 
a plateau much earlier than before, immediately after to (which is also set to 4). Of course, the excited state(s) is(are) 
much noisier and we could not reliably extract information from that, but if one is only interested in the ground state, 
this is not an issue. However, one could apply the GPoF method to a large correlator basis to perhaps get cleaner 
signals for the lower lying states. On the right we perform the shift to the two operators, creating a 4 x 4 basis, but 
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FIGURE 3. This is the same as in Fig. 2, but with a single shift of T = 4 for the operator(s) included. 



after the SVD cut, this is reduced to three operators. In either case, we see that the ground state is unchanged, and a 
plateau is reached far sooner than without the use of the GPoF method. 

For the unshifted operator if we perform a correlated fit to the correlator using an exponential in the range 
t G {29,40}, we get a t m^ = 0.2770(50) with a confidence level of 78%. For the single shift of z = 4, we use the 
larger (and earlier) range t s {5,28}, and get a t m^ = 0.2825(17) with a confidence level of 84%. Adding more shifts 
does not change this result, but it is amazing how well the signal improves with just this single shift. 

It can be seen what is happening qualitatively, since the correlator C(t + x) comes from the correlation with an 
operator that is T time steps away from the other operator, and thus excited states are not going to contributed much to 
that, and even less to C{t + 2t). The ground state, however, is going to still contributed significantly to this correlator, 
so its signal, in a sense, is effectively enhanced. 

In the three-point case, working through the correlator matrix, we would get 

<?-*{ti,t,tf) C^(t U t,tf + T) \ 

C 3 -^(ti,t + T,tf + T) C 3 -P t (f,-,f + T,f / + 2T) J ■ K ' 

So unlike the two-point function, we have to actually use three different sink locations, and thus generate a factor of 
three in the computational cost. However, if this allows a better determination of our ground state signal, this should 
be worth the additional propagator generation. 

Once the two-point correlator matrix is diagonalized with the vectors x, we form a matrix Vy = (x ; ), and then use 
that to diagonalize the three-point correlator 

C^%t,tf) = V- l C 3 -v\t„t,t f )V . (12) 

Currently, the analysis using this approach is under way and we hope to see a noticeable improvement as we did in the 
two-point case. 
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